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1. Introduction 

The interest in the impact of the cosmological expansion on bound astrophysical system 
can be traced back to the paper by Einstein and Straus [Ij. Over the subsequent years 
different approaches have been used to gain some insight into the interplay between 
cosmology and astrophysical structures [21 EH HI El E] . With the advent of the discovery of 
the accelerated Universe [7] , it seems timely to pick up the topic once again and contrast 
the effects of different models which can explain the current accelerated phase. Indeed, 
in [8] such an examination was performed for phantom [9] and quintessence cosmologies 
[TO] . In the present article we will use the same approach as [2] and [Sj and apply 
it to models with a positive cosmological constant and Chaplygin gas models. These 
two models will serve us as representatives for cosmological models which can explain 
the current stage of acceleration. In addition we will also invoke a general model with 
variable equation of state of the form uj x = p x / p x (DECDM). The approach mentioned 
above consists in generalizing the Newtonian limit by allowing terms which are specific 
to the background cosmology. This is not limited to the inclusion of the cosmological 
constant. Notice that the effects of the background cosmology on bound systems are now 
in general time dependent. This is, however not in contradiction with the fact that we 
rightly consider most of the astrophysical structures to be in gravitational equilibrium. 
The time dependent background cosmology will not tear the system apart, unless its 
density is diluted, but change only, over cosmological time, the internal properties of 
the system (like inner velocities or angular velocity etc). As a mathematical tool we will 
make use of the virial theorem ([TTJ for different methods) by allowing certain internal 
properties to be time dependent. The virial theorem can be easily derived even in the 
presence of the cosmological background as the latter enters only the Newtonian limit 
(i.e. the Poisson equation for the gravitational potential which includes now a time 
dependent part). 

As shown in JT2] Dark Energy can affect certain astrophysical static properties 
which have to do with either the virialization of the system [13J or the motion of 
test bodies GO. For instance, the cosmological constant sets itself scales of distance 
(rA = A -1 / 2 ) time and mass, which are of the same order of magnitude as the radius 
of the observed universe, the age of the universe and the total mass of the universe (the 
so called coincidence problem). At the first glance it looks hopeless to expect any effects 
from A at astrophysical scales. This is so indeed as long as no other scale enters the 
theory. If the latter appears in the theory a combination of the large cosmological and 
the small internal scale can yield values of astrophysical relevance (e.g. (t^a) 1 / 3 where 
r s is the Schwarzschild radius) [Hj. Furthermore, general conditions for detecting the 
effects of variations of 'constants' represented by scalar fields (just like the cosmological 
constant and quintessence models) at scales ranging from terrestial scales have been well 
established in [T5l [T6l [T7] . Another enhancement mechanism of background cosmology 
emerges when we are dealing with non-spherically symmetric objects. In such case 



| We work in the so-called geometrized units, i.e. Gn = c = 1 
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the ratio of two length scales to some positive power (bigger than one) goes often 
hand in hand with A. For this reason, in examining time dependent effects, we will 
mostly employ ellipsoidal configurations. For general consideration of equilibrium in 
the spherical case see [18] and [19], the quasi spherical case has been discussed in [20] . 

In this paper, we will entirely concentrate on effects of Dark Energy followed over 
cosmological times, past and future. We do not expect these effects to be large, however, 
we think that it is of some importance to probe into these matter, especially in view of 
the accelerated universe. In [21] we restricted ourselves to effects over small time scales. 
In this sense these two articles complement each other. 

2. Local dynamics with background cosmology 

In standard cosmology, the universe is described as an ideal fluid which implies the 
properties of homogeneity and isotropy. Such properties are mathematically represented 
in Friedman- Robertson- Walker line element, which is written for a flat universe as ds 2 = 
— dt 2 + R 2 (t) (dr 2 + T 2 dQ 2 ). Einstein field equations and the FRW line element yields 
the acceleration equation and the well known Friedman equation, written respectively 
as 

= H{tf = | wpb (t), |j| = --7T [p h (t) + 3 Ph (t)\ , (1) 

In conventional cosmology, the total energy density pb is a contribution of a radiation, 
cold dark matter and dark energy. By neglecting interaction among these components, 
each one evolves through the conservation equation p = —3H(p + p), so that radiation 
scales as p ra( j = p Ta ,d(R = 1)-R -4 , cold dark matter scales as p cc i m = Pcdm(R = ^)R~ 3 , 
while the Dark Energy component associated with the cosmological constant is p vac = 
constant. The major contribution in connection to the cosmological constant is 
approximately 70% of total energy density [TJ. 

2.1. Dark energy 

As pointed before, almost the 70% of the content of energy density in the unverse is 
ruled by a dark energy component, and astronomical observations imply that at the 
present time this dark energy component might be represented by the cosmological 
constant. That is, at the present time, the equation of state that displays the best 
fit with astronomical observation is simply p vac = —p vac , i.e, u) — — 1. Nevertheless, 
several models has been proposed in order to reproduce the present value of dark 
energy but maintaining a time dependent energy density. One generalizes the equation 
of state for the dark energy by p x [R) = u x (R)p x (R). Dark energy then scales as 
p x (R) = p x (R = l)a~f( R >, where the function f(R) is defined as 



R(t) 

R(t) 
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The model is constrained by requiring that at the present time the dark energy density 
reproduces the vacuum energy density associated with the cosmological constant, that is, 
p x (R = 1) = p vac . As pointed before, the case u x = — 1 corresponds to the cosmological 
constant p x = p vac = A/87T, while for > u x > — 1 one stands in the Dark Energy 
realms. Models with uj x < — 1 lead to future singularities in the so called Phantom regime 
[9]. The description of a dynamical dark energy is compatible with the description of 
inflationary models. 

In the DECDM model, the Friedmann equation and the acceleration equation are 
written as 

= H o Q ^MR) ^ = -Hfa^h 2 (R) = ~irp^h 2 (R).(3) 
The functions h\ }2 (R) are given as 

/il ( i2 ) = ^ ijR -3 + jR -/(i?) ) (4) 



'vac 



h 2 (R) = \ 



R- 3 + R- f{z) (1 + Su x (R)) 



^cdm 
^vac 



2.2. Chaplygin Gas 

The main stages the universe has passed through in the standard cosmology scenario can 
be reproduced by the introduction of the equation of state of an exotic ideal relativistic 
gas written as 

Pch = ^lPch tT, (5) 

Pch 

The relevance of this equation of state lies in three different facts: a) it violates the 
strong energy condition, which is necessary to obtain an accelerated phase at the present 
time b) it generates a well definite speed of sound, which is relevant for the process of 
structure formation and c) it unifies the early radiation or dark-matter behavior with 
the late dark energy dominance. The so-called pure Chaplygin gas is obtained with 
Ki = and 7 = 1. 

The time evolution of the Chaplygin Gas-energy density is given from the integration 
of mass-energy conservation equation as 



Pch (R)= [A + BR-^l", (6) 
with n = 3(1 + Ki), (3 = 7 + 1 and where A and B are integration constants given as 

B = Pch (R = l) 13 - A, A = -k 2 = -a^l 13 , (7) 

n n 

with p c h(-R — 1) — Pcrit as the Chaplygin Gas energy density at the present time (only 
for a flat universe). In order to write the integration constants as we have done, we 
have re-parametrized the Chaplygin Gas equation of state with k 2 = ct n hl]f such that a 
finite age and a present accelerated phase for the Chaplygin gas universe is reached for 

n { 3\ p 

"■ < s(sf) • (8) 
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This bound can be derived after integration of the Friedmann equation 



H(R) 2 = % = ^ 



A + BR 



-n/3 



(9) 



R 2 3 

in order to determine the age of the universe. The age of the universe is given in terms 
of hypergeometric function as 

3 \P 3 [T 
2^ 



T 



qH 
x 2 F l 



57T 



-Ct r , 

n 



(10) 



q 1 2n(3 + n + A 3 
2n/3' 2/3' 2^/3 ' ~n Qr 



3\P 3 
— --«, 

57T/ n 



with q = 3u + 7. For the pure Chaplygin gas one then obtains a finite age and the 
observed acceleration at the present time only if the parameter a is such that < 10~ 2 , 
which in turn implies a bound on the Chaplygin parameter k 2 < pl rit . This warranties 
that the equation of state is a; c h(today) > — 1. 

The solution given in Eq.Q allows us to describe the behavior of Chaplygin gas at 
different ages. For early times the Ki-term dominates and Chaplygin Gas behaves as 
p c h = B l IPR~ n . This scales as radiation for n = 4 (ki = 1/3). Hence we could 
in principle relate the integration constant B with cosmological parameters of the 
concordance model as B = pf ad (-R = 1). For K\ — 0, the model displays a matter 
dominated epoch at early times with p c h w B 1 / /3 R~ 3 . In this case one can write the 
integration constant B = Pmat(-R = 1) an d also the constant A as A = pf rit (l — ^mat)- 
Finally for large times, Chaplygin gas acquires a vacuum-like behavior with 

"3a n ^ 



Pchif 



oo 



A l/P 



n 



Hi 



(11) 



This would lead to an effective vacuum energy density p^ c which can be associated to 
the current vacuum energy density p vac via 



r vac 



Pv 



i + 



B 
A 



(12) 



With pl^ c > p vac , the effects of this effective vacuum energy density may be relevant in 
the Newtonian limit when we explore the equilibrium conditions of large scale structures. 
Let us concentrate on the case K\ = 0, which does not reproduce an early radiation 
dominated era but still unifies an early dark matter dominance and a vacuum dominated 
era at large scale factors. The acceleration equation for this model reads as 



W) = ~ t nPch{R)rich{R) 



(13) 



4 

-71 

3 



A + BR(t) 



-3/3 



1 - 



3k 2 



A + BR(t) 



-3/3 



where 7/ c h = 1 + 3c<j c h and 

^ch = ^ch(-R) = 



Pch 
Pch 



-K 2 p c h{R)~ 



(14) 
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Figure 1. Left: Equation of state for Chaplygin gas w c h = Pch/pch for different 
Chaplygin parameters. Right: Acceleration equation as a function of the (future) 
scale factor for different Chaplygin parameters. 

is the time dependent equation of state. In Fig. [T]we show the behavior of the equation 
of state u c h and the acceleration equation for different values of Chaplygin parameters. 



2.3. Newtonian description 

The effects of the background can be explored through the Newtonian limit of field 
equations [22], from which one can derive a modified Poisson's equation [2] 



V 2 $ = 47r<5p-3 - , (15) 




where Sp is the overdensity that gives rise to the "gravitational contribution" of the 
potential. The total energy density within the clustered configuration p is a contribution 
of the background pb and the collapsed fraction Sp. The background energy density 
depends on the model we are interested in: 

(m A J^dmOR) DECDM 
p b = p(R) -Sp=t . (16) 

{p c h{R> Chaplygin 

Note that the last expression implies that in the DECDM model one has assumed that 
only the cold dark matter component is what collapses and forms vitalized structures, 
while in Chaplygin gas the perturbation is done on the total energy density. The solution 
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for the potential can be simplified as 



*(r, R) = $ grav (r) - ( - ) r\ $ grav (r) = - / -^ydY, (17) 



A 

where one has neglected the terms associated to boundary conditions (the potential 
cannot be zero at infinity) |22j. The last term of Eq.( |T5|) reduces to p vac for u x = — 1 and 
negligible contribution from the cold dark matter component, that is, for low red-shifts. 
This limit represents the Newton-Hooke space time, which aside from the gravitational 
interaction has an external force of the Hooke form ~ Ar 2 . 

Gravitational equilibrium is represented through Euler equation pii = —Vp — pV^, 
where p is given by Eq.( jT6l) . The next generation of equilibrium equations comes from 
taking moments on Euler equation and derive the second order virial equation 

<l " Ji!: 2K ih + Wg &w 4 311,, -i I ■- I J,,. (1> 




dt 2 ^ k ^" y iK 

where Wf^ av is the gravitational potential energy tensor, whose trace corresponds to 
the gravitational potential energy, K.^ is the kinetic energy tensor, X^ is the moment 
of inertia tensor and IL^ is the dispersion tensor. The trace of these quantities lead to 
the well known forms for gravitational potential energy: 



W srav = -] pr^$ grav d d r = -J p$ grav d d r, (19) 

together with X = J v pr 2 d 3 r and II = J v pd 3 r. As usual, the set of equilibrium equations 
is closed with the equation for mass conservation, energy conservation and an equation 
of state p = p(5p, s). 

The tensor virial equation can be thought as a differential equation for the moment 
of inertia, which in turns is converted to a differential equation for the parameters 
determining the geometrical properties of the configuration. In order to explore the 
contribution from the background, let us write the virial equation (fIBl) as 

d 2 X ik 2 




(20) 



dR 2 R 2 H 2 (R) 
3. Dynamical equilibrium 

If we assume equilibrium via X « 0, we obtain the virial theorem 

| W grav| = 2/c+ f^j X) IC = T + hl. (21) 

Strictly speaking, a formal equilibrium configuration is never reached since the energy- 
like terms in <HM are time dependent. Nevertheless, we can use equation (|2"TI) by 
assuming that as long as R evolves in time, the configuration evolves through successive 
states of equilibrium. 

The expressions derived in the last section, specially ( |T5l) and ( l21i) can be used for 
testing dark energy models on configurations in equilibrium. In this section, we will 
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concentrate on possible effects of the two cosmological models described above on 
cosmological structures with non-spherical symmetry, such as low density galaxies and 
galactic clusters. 

For ellipsoidal homogeneous configurations, the gravitational potential tensor and the 

moment of inertia tensor are written as 

s 

:7r 2 p(p - p h )aia 2 a 3 a^A i 8 ik = -2rt:5pA i l 1 



ik 



15 



ik ■ 



In 



4 



--'Kpa l a 2 a 3 a 2 5 ik , 
15 



(22) 



where the quantities A4 are functions of the eccentricities of each case: for an oblate 
configuration we have a\ = a 2 > a 3 , el hl3ute = 1 — q 2 , while for prolate a± = a 2 < a 3 , 
e proiate = 1 — <h 2 '' ■> with qi = di/ a\. The functions Aj are given as [24J: 



A 1 = A 2 



arcsm e 



1 - e 2 r 1 1 , 

n m 

1-e 2 2e 



1 + e 



Oblate 
Prolate 



(23) 



and 



A, 



VT 



2e VI -e 



arcsm e 

e 

- 1 



Oblate 



Prolate 



(24) 



In the limit e — > one has Ai 2/3 for the prolate and oblate configurations. Using 
Eqs.f l22p . the potential energies are then written for the models of interest as 



W ik = ^grav + (25) 

ri 

8 4 / 2 \ 

-—-K 2 p 2 aia 2 a 3 a 2 i A i 5i k - -nr] x p x l ik + 2np cdm \ Ai - -J l ik 

8 / 2 \ 

-— ir 2 p 2 a 1 a 2 a 3 a 2 Ai5 ik + 2rcp ch (a { - -r] ch ) l ik 

for DECDM and Chaplygin Gas respectuvely, with r/ x = 1 + 3c<j x and p x = p vac (l + z)^). 
Note that the last term on the DECDM model vanishes in two different situations: 
the first is in the case of spherical symmetry, and when we neglect the contribution 
from the cold dark matter component with respect to the proper density of the system 
and the dark energy contribution. This approximation is only valid for very low red- 
shifts. At high red shifts we approach a matter dominated universe and hence the term 
proportional to p c d m is relevant. 

3.1. Oblate systems 

When considering oblate configurations {a\ = a 2 > a 3 ), we assume that the kinetic 
energy is due to constant rotation along the minor axis. We then may write the kinetic 
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energy tensor as 

Tik = ^ {^vot-^ik ^roti-^fcj'^rotj 



(26) 



In order to determine the angular velocity, one may consider a non-zero isotropic 
dispersion tensor = 5^11. We then have 



n 



I Ws rav I 

I ' * XX 



rav 

zz 



rot 



T 

■J-xx 



R 



1 



T 

-'-zz 

J~ 

-J-'.r.r 



(27) 



which follows from eliminating the trace of the dispersion tensor from the tensor virial 
equations. Using Eq. (l2"5|) . we can explicitly write for the angular velocity 



n 2 



2n5p 



A 1 -A 3 qi) x 



(28) 



/ 



V 



2n5p(A 1 - A 3 q 2 3 ) 



-f^0^vac^2(-R) 

x I 4 

-npch(R)vch(R) 



DECDM 
Chaplygin 



In the Newton-Hooke space time with u x = — 1 (h 2 — > — 1) this reduces to 



(29) 



which in turn reproduces the Maclaurin formula for p vac = 0. 



3.2. Prolate systems 

For prolate configurations (a± = a 2 < a 3 ), we solve for the velocity dispersion of main 
components rather than for an angular velocity (rotation of prolate configurations, 
although rare, has been observed and some properties of this situation has been explored 
in [2T]). Hence we use the virial theorem (]2"Tj) with a the kinetic energy given by 

\j v p{v 2 )dh = ^K P 4^ 2 ). 

In analogy with Eq. (|28|) . the virial theorem then implies for the velocity dispersion 

/ 



K 



(30) 



(v 2 ) 

lllbp 



~a\ {2A X + q\A 3 



1 + 



2 + qi 



\ 



27r8p(2A 1 + q 2 A 3 



HlVt Y3iC h2(R) 

4 | (31) 

-7rp ch ( J R)^ch( J R) 



for the DECDM and Chaplygin Gas model respectively. As in (128]) . one obtains a scale 
factor (or red-shift) dependence of the velocity dispersion through the contribution of 
the cosmological background. In the Newton-Hooke space-time one reduces to 



(v 2 ) 



-a\{2A 1 + q 2 3 A 3 



4 Pvac 



2 + qi 



3 Sp \2A 1 + q 2 3 A 3 



2n5p 5 

which in turn reduces for spherical symmetry to the known expression 
(v 2 ) 1 



2n5p 



m —n 
5 



1 - 2 



Pvac 
lp~ 



(32) 



(33) 



From Global to Local Dynamics: Effects of the Expansion on Astrophysical Structures!^) 




0.985 



0.6 



1.6 

Scale factor 



2.6 



0.9985 



0.6 



1.6 

Scale factor 



2.6 



Figure 2. Equation (|34|) for A„, for different densities and eccentricities in a ACDM 
model. 



The structure of the resulting expressions for the angular velocity and the velocity 
dispersion can be written in a similar way. Let A measures the ratio of each velocities 
with respect to the velocities when there is not background contribution. We then may 
write 



A 



rot 



1 



(34) 



I (v 2 )( Pb = 0) 
where the geometrical factor Q is written for each 



2 + g| 
2A 1 + qlA 3 



Oblate, rotational velocity 



Prolate, velocity dispersion 



(35) 



Returning to the oblate case, in a generalized Newton-Hooke space-time (i.e, with 
— 1 < uj < ), one observes that the Dark Energy component decreases the angular 
velocity (or the velocity dispersion in the prolate case) in order for the system to 
maintain equilibrium (since 1 + 3uj x < 0). At the present time (R = 1, z = 0) we 
get A n ~ 1 - 0.8#(e)(p vac /5p), where 



9{e) 



-e 
3 



(3 - 2e 2 



arcsm e 



3e(l-e' 
0, which implies Q 2 



(36) 



An extreme situation is reached in the limit A^ — > 0, which implies Vt z — > 0. For e ~ 0.8 
this would require ( = 2p vac /<5p ~ 0.5, which is a very diluted configuration. For higher 
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Figure 3. Equation ([34| for A„, for different densities and eccentricities in a Chaplygin 
Gas model for typical values of Chaplygin parameters. 

eccentricities (say e ~ 0.97) one would need a system with ( ~ 0.2 which is still a 
rather low value. In the ACDM model Q = is reached at a red-shift given by [26] 
z c = [(2f2 vac ((g(e) — 2))/(Cg(e)f2 c d m )] 1 ^ 3 — 1- For realistic examples as 5p ~ 200p cdm the 
value z c could be reached for e very close to 1. That is, the value z c requires extreme 
flat objects. 

In figures [2] and [3] we show the behavior of the ratio for different values of eccentricities 
and proper densities in the DECDM case with u x = — 1 and in the Chaplygin Gas 
model, respectively. The larger effects appears for low densities and large eccentricities, 
as expected. However, the fractional change is very low. For 5p = 10 2 p vac , the larger 
effect occurs for scale factors > 2.6 with a change of 1%. 

In Fig. H] we show the behavior of Ar> — 1 as a function of the red-shift for different 
densities and e = 0.9, 0.95 for three different values of the equation of state for Dark 
Energy. It is clear that the effects associated with a cosmological constant are stronger 
than the ones associated with other Dark Energy models, but in general those effects 
are small for realistic values of ( as used in the figure. We would have to measure the 
angular velocity at different red-shifts very exactly to see an effect over a range of z. 
However, the difference between the models is more significant. 
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Figure 4. Equation (|34|) for A si — 1 in terms of the red-shift (up to a ~ 0.7), for 
different densities and eccentricities in a DECDM model for typical values of equation 
of state. [26] 



4. Dynamical evolution 



In this section we want to explore the virial equation with the cosmological background 
without insisting on the equilibrium condition (dynamical equilibrium). Even though 
we will not make extensive use of these new equations, we want to demonstrate that the 
equilibrium condition is not a necessary ingredient while studying the effects of different 
background cosmologies. The results are differential equation for the geometrical 
parameters of a given configuration. Our starting point is equation (|20l) . In the first 
approach, we will concentrate on pressurless systems with homogeneous density. The 
problem is the reduced to determine ax = ai(-R) and a 3 = a 3 (R). As pointed out before, 
Eq. ( l20i) represents a set of two non-linear coupled second order equations for i = x,z 



d 2 Iu 



dR 2 R 2 H 2 (R) 



2T U + wr + (I 




ldX, 



2 dt 



(37) 



where we have used Eq. fl25j) . The derivatives of the moment of inertia tensor are given 
as follows 



d2u 
dR 

d 2 2ii 
dR 2 



8 2 dai ^da; 



15 



dai \ da; 
a,- — + a,- — 
v dR dR 



2aJ 



da; N 
'dR 



+ a 



d 2 &i 
dR 2 



(3f 



X 
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These expressions enclose our approximation: the volume of the system is constant even 
if the semi- axes may change, so that one can write the moment of inertia tensor as 

X ik = ^Ma 2 6 ik , (39) 
where the mass M is 

4 



M 



M b + 5M = -7rala 3 p h (R) + 5M } 8M = f 5pd 3 r. (40) 
3 Jv 

In order to deal with the kinetic energy tensor, we follow [27J and transform from 

proper coordinates to fixed co-moving coordinates via = A ia x a , where A ia is such 

that Ai a = ai5i a and hence the ellipsoid is characterized by the constraint b a px a x^ = 1. 

We have 

1 r fdA ia \ /dA k/3 x 



T * = 2'J r [-ir){-£ L )*"> A **< (41) 

where A = det(A a p) = aia 2 a 3 . Using the definition of the moment of inertia tensor one 
has 

X ik = pA ia A k/3 A J x a xpd 3 x , (42) 

and hence 

T, = \ A Zt A te ( d 4^)( d -¥) **■ (43) 



We then may write 



R . k($ )V.<*(S0V.j|^ m 

for the DECDM and Chaplygin Gas respectively. Equation fl37j) together with fl38l) are 
the differential equation to be solved. The functions Ai depend on the eccentricity of 
the system and hence they are also function of the semiaxis Ai(z) = Ai(ai(R) , a 3 (R)) . 
The differential equation can then be cast into the following form 

I • (45) 



d 2 ai 



dR 2 R 2 H 2 (R) 



2T U + Wu fR\ 1 da ; 



la \R ) Oi dR 



daj 
dR 



5. Conclusion 



In this article we have examined the effects of the expansion of the Universe on certain 
quasi-static properties of large astrophysical bodies. We have done this by invoking a 
dynamical equilibrium i.e. the time dependent response of the astrophysical objects to 
the time dependent cosmological background is taken into account by allowing some 
internal properties to become epoch dependent. In spite of the fact that the expansion 
of the Universe is accelerated the effects are rather small, at the most few per cents over 
a large cosmological time stretch in the case of the Chaplygin gas model. However, a 
qualitative investigation, done in the present paper, seems mandatory to establish the 



From Global to Local Dynamics: Effects of the Expansion on Astrophysical StructureslA 

size of the effects. Different equations of state can be clearly distinguished theoretically, 
however, in practice this would require a good knowledge of the angular velocities at 
different redshifts. In spite of the smallness of the effects, the results clearly demonstrate 
that background cosmology affects, in principle, local properties of astrophysical bodies 
and furthermore, the details of these effects are model-dependent. In this paper, we 
focused on the the dynamical equilibrium. One can enlarge the concept of the influence 
of background cosmology by dropping the equilibrium condition and go over to a fully 
dynamical scenario. A first step in this direction was done on section 4. Supplementing 
these equation with Euler equation for the angular velocity, one has a full set of 
differential equations determining the time evolution of the parameters during (and 
after) virialization. We will come back to this point in future publications. 

References 

[1] Einstein, A. Straus E. G., Rev. Mod. Phys. 17, 2 and 3, 1945 

[2] Noerdlinger P., Petrosian V., Astrophys. J, 168, 1 (1971) 

[3] Bona C, Stela J., Phys. Rev. D36, 2915 (1987) 

[4] Baker G., A., arXive |astro-ph/0112320| 

[5] MacVittie G., C, Mon.Not.Roy.Astron.Soc. 93, 325 (1933) 

[6] Nunes, N. J. and Mota, D. F., Mon. Not. Roy. Astron. Soc. 368:2 751 (2006), |astro-ph/0409481 
[7] Perlmutter S. et al, Nature 391, 51, (1998); Perlmutter S. et al., ApJ, 517, 565, (1999) 
[8] Nesseris S., Perivolaropoulos L., Phys. Rev. D70, 123129 (2004) 

[9] Nojiri S., Odinstov S. D., arXiv: |hep-th/0505215| Godlowski W. , Szydlowski M., arXiv: 

|astro-ph/0507322| Polarski D., Ranquet A., arXiv: |astro-ph/0507299[ Hu W., Phys. Rev. D71, 

047301 (2005). arXiv: |astro-ph/0410680| 
[10] Barreiro, T., Copeland, E., Nunes, N. J. Phys. Rev. D 61, 127301 (2000); Steinhardt, P. J., Wang, 

L., Zlatev, I., Phys. Rev. D59, 123504 (1999); Albrecht, L.R., Skordis, K., Phys. Rev. Lett. 

84 2076 (2000); L. Wang, R. R. Caldwell, J. P. Ostriker, Paul J. Steinhardt. ApJ. 530, 17-35 

(2000), arXiv: |astro-ph/9901388| 
[11] Cardoso, V. and Gaultieri, L., Class. Quantum Grav. 23, 7198 (2006) 
[12] Balaguera-Antolmez, A., Nowakowski, M., Astron. Astrophys 441, 23 (2005) 
[13] Mota, D. F. and van de Bruck, C, Astron. Astrophys. 421, 71 (2004) 

[14] Balaguera-Antolmez, A., Bohmer, C, Nowakowski, M, Class. Quantum Grav. 23, 485-496 (2006) 

[15] Douglas J. Shaw, John D. Barrow Phys. Lett. B639 596-599 (2006) 

[16] Douglas J. Shaw, John D. Barrow Phys. Lett. B73 123506 (2006) 

[17] Douglas J. Shaw, John D. Barrow Phys. Lett. B73 123505 (2006) 

[18] Boehmer, C. C, Gen. Rel. Grav 36, 1039 (2004) 

[19] Boehmer, C. G. and Harko, T., Phys. Rev. D71, 084026 (2005) 

[20] Debnath, IT., Nath, S. and Chakraborty, S., Mon. Not. Astron. Soc. 369, 1961 (2006) 
[21] Balaguera-Antolmez, A., Mota, D., F., Nowakowski, M., Class. Quantum Grav. 23, 4497-4510 
(2006) 

[22] Nowakowski, M., Int. J. Mod. Phys., D10, 649 (2001) 

[23] Nowakowski, M., Sanabria, J.-C, and Garcia, A., Phys. Rev. D66, 023003 (2002) 
[24] Binney, J & Tremaine, S., Galactic Dynamics, Princeton University Press, 1987 
[25] Wang, L., and Steinhardt, P. J., ApJ 508, 483 (1998) 
[26] Balaguera-Antolmez, A., Nowakowski, M., arXive: astro-ph/0603624 

[27] Peebles P.J.E., The Large-Scale structure of the Universe, Princeton University Press, 1980. 



